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Abstract 

As a preparation for the numerical study of the SU(2) gauge theory with gluinos, 
the spectral properties of the fermion matrix and the masses of pseudoscalar and 
scalar states are investigated in the quenched approximation. The behaviour of the 
disconnected fermion diagrams on small lattices is also studied. 



1 Introduction 

The most remarkable developments in non-perturbative Quantum Field Theory (QFT) in recent 
years have been triggered by the exact solution of certain supersymmetric gauge theories by 
Seiberg and Witten [p]]. These exact results are exploiting the electric-magnetic duality and 
imply rigorous proofs of confinement and chiral symmetry breaking. They open the door for 
a new non-perturbative procedure: after solving the theory in a highly symmetric point of 
the parameter space one can calculate the renormalized effective action by deformations or by 
expansions in symmetry breaking parameters. (For examples, see ref. [0, [|.) 

The conventional approach of non-perturbative QFT based on lattice regularization is in 
some sense opposite. There, in the cut-off theory, the symmetries are broken, with the important 
exception of exact local gauge symmetry, and have to be restored in the continuum limit by 
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parameter tuning. In spite of the great successes of the Seiberg-Witten approach, it would be 
important to establish a connection to lattice regularization by investigating supersymmetric 
gauge theories on the lattice. A possible way to do this was put forward some time ago by 
Curci and Veneziano ||J]. (For the case of N = 2 supersymmetric gauge theories considered in 
see also §.) 

The simplest supersymmetric gauge theory is pure SU(2) Yang-Mills theory with massless 
gluinos. In case of massive gluinos the non-vanishing gluino mass softly breaks supersymmetry. 
Recently some possible fermion simulation algorithms for theories with gluinos have been devel- 
oped and tested || [7]. A numerical simulation with dynamical gluinos is in progress || |J. As 
a first preparatory step, so called "quenched" studies are useful, where the gauge configurations 
are prepared in pure gauge theory. 

In the present paper we investigate some interesting features of the SU(2) gauge theory 
with gluinos in the quenched approximation. The spectral properties of the fermion operator 
are studied, which are important in the class of local bosonic fermion algorithms [Tt| we shall 



apply. First information about the interesting region of the fermion hopping parameter, where 
hadronic masses are small in lattice units, can be obtained from determining the masses of 
pseudoscalar and scalar bound states. In our case these are pion-like, sigma-meson-like or 
eta-meson-like, states which are held together by the colour forces among the gluinos in the 
adjoint representation. (In case of SU(2) this means colour triplets.) Therefore, we shall call 
them adjoint-pion (shortly a-pion), adjoint- sigma (shortly a-sigma) and adjoint- eta (shortly 
a-eta), respectively. For the a-eta state, the computation of the disconnected fermion diagrams 
is necessary. 

In the next section we report on our quenched calculations carried out at the bare gauge 
coupling parameter (3 = A/g 2 = 2.3 , which roughly corresponds to the lower edge of the scaling 
region in pure SU(2) lattice gauge theory with Wilson action. 

2 Numerical results 

2.1 Spectral boundaries of the fermion matrix 

The fermion matrix Q relevant for gluinos can be defined as 

4 

Qyv,xu — Qyv,xu[U] — Syx^vu K ^ ] ^?/,a;+/i(l T/i)^u«,x^t &y+fi,,xiX Tm) ~^vu,y[i " (-Q 

Here, as usual, x,y,... denote lattice sites, u,v,... colour triplet indices, fi the unit vector 
in direction fi = 1,2,3,4 and K the hopping parameter determining the gluino mass. The 
gauge link variable in the adjoint representation V X41 is connected to the usual fundamental link 
variable U x ^ by the relation 

v rs,xn = V rs , x ,[U] = 2Tr(Ul%U xp T s ) = V r * s ^ = V~Z , (2) 
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if T r is a group generator (for SU(2) T r = |r r is proportional to the Pauli matrix r r ). 

In local bosonic updating algorithms the spectral properties of Q and Q^Q are important, 
because the polynomial approximations defining the bosonic action are optimized on the eigen- 
values. In the two-step local bosonic algorithm |J the polynomials are optimized in an interval 
[e, A] containing the eigenvalues p of Q^Q. For a hermitean matrix as Q^Q the extremal eigen- 
values p min = minp and p max = maxp can be obtained by known iterative methods [|11]] (for 
recent work on the spectrum of the Wilson-Dirac operator in QCD see [0). In the same way, 
the extremal eigenvalues of the hermitean matrices 

Qr = 7^(Q + Q ] ) , Q^i(g-gt) ( 3 ) 

can also be determined. Since for vectors v normalized to unity we have 

Q = Qr + iQi , Re (v\Q\v) = (v\Q R \v) , Im (v\Q\v) = HQ» , (4) 

the spectral boundaries of Q itself can also be estimated. For instance, if A denotes eigenvalues 
of Q and Xr eigenvalues of Qr, we have 

\Rmin = minAij = minRe ) < minReA , (5) 

and similarly for the maximum \R max of Xr and the extrema \i m in,max of A/. The numerical 
simulation results for the quenched averages of these extremal eigenvalues at (5 = 2.3 for 
0.180 < K < 0.205 on 8 3 ■ 16 lattices are shown in table |l|. The distributions of p m i n and 
^Rmin for K = 0.200 are also shown in figures |T] and |^, respectively. The statistics is based 
on at least 320 independent gauge configurations per hopping parameter value. The iteration 
for the extremal eigenvalues was stopped if the relative deviation of the new estimate from the 
previous one was less than 10 -3 . 



Table 1: The quenched averages of the extremal eigenvalues defined in the 
text on 8 3 -16 lattice at (3 = 2.3 . The statistical errors are given in parentheses. 



K 


Pmin 


Pmax 


^Rmin 


^Rmax 


A/max ^Imin 


0.180 


0.00879(2) 


4.4873(6) 


-0.0801(1) 


2.0690(2) 


0.86357(5) 


0.190 


0.00434(2) 


4.764(1) 


-0.1375(2) 


2.1287(3) 


0.91151(9) 


0.200 


0.00197(1) 


5.050(1) 


-0.1958(2) 


2.1885(3) 


0.95948(9) 


0.205 


0.00138(2) 


5.194(1) 


-0.2245(2) 


2.2178(3) 


0.98373(9) 



As table [I] and figure |1] show, the lower limit of the spectrum of Q^Q is of the order of 
10 -3 , and the ratio of the largest to smallest eigenvalues is 10 3 -10 4 . This is the range one also 
encounters in "unquenched" numerical simulations with dynamical gluinos [[J . It is interesting 
to observe that the lower limit of the spectrum of Qr defined in (|3|) is negative in the whole 
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range of the hopping parameter considered. Although the minimal real part of the eigenvalues 
of Q might be somewhat larger, it is probably also negative. This shows that in the relevant 
hopping parameter range the eigenvalues of Q do not remain in the half plane with positive 
real part, but surround zero in the complex plane from all sides. 

2.2 Masses of a-pion and a-sigma 

The connected contributions to pseudoscalar and scalar bosons are 

Tr {l.QyHQ^y} , Tr {Q- 1 ^- 1 } , (6) 

respectively. These give the correlation functions of the corresponding flavour non-singlet 
mesons in models with at least two (adjoint) fermion flavours (N a > 2). Such states, which we 
call a-pion and a-sigma, respectively, do not exist in the gluino model, which is equivalent to 
an (adjoint) flavour number of N a — -z. Nevertheless, the masses of the a-pion and a-sigma are 
interesting, because the hopping parameter region where they become small in lattice units is 
the interesting one, where also other physical colour singlet hadrons become light. 

Similarly to the pion in quenched QCD, one expects that the a-pion mass-squared vanishes 
linearly as a function of 1/K. This is the remnant of chiral symmetry for N a > 2 in the 
quenched approximation, because the bare quark mass in lattice units is m q = \{1/K — 1/Kq). 
Here K = K (j3) is the hopping parameter where the pion mass vanishes. Of course, since the 
gluino corresponds to N a — ^, in the gluino model there is no analogous chiral symmetry and 
the exact value of the "critical" hopping parameter (K cr ) where the a-pion has zero mass is not 
directly relevant. It is interesting just as an indicator of the scaling region where the continuum 
limit has to be performed. 

The a-pion mass at (3 = 2.3 as a function of 1/K is shown in figure |3]. The linear extrap- 
olation works very well and gives K cr = 0.2151(3) . The masses of a-pion and a-sigma are 
displayed together in figure |j. As expected, the a-sigma mass is higher and does not seem to 
vanish at K cr . The results in figures |^ and [| were obtained on an 8 3 • 16 lattice. The masses 
were always extracted from timeslice correlations by a fit 

c + c 1 {e- mt + e- m ^} , (7) 

with T being the lattice time extension. For comparison, in figure [| at K = 0.180 also the 
masses on 4 3 • 8 are included, which do not deviate too much from the corresponding masses 
on 8 3 ■ 16. This shows that, not too close to K cr , even a 4 3 ■ 8 lattice can give first indications 
on the mass spectrum. 
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2.3 Disconnected contributions: the mass of a-eta 



The correlation function of the pseudoscalar meson made out of Majorana gluinos can be 
obtained from the general formulas given in ref. |J: 

^7fe*»)(^-7fe*.)> = (Tr {j 5 Qyy}Tr { l5 Q^} - 2 Tr {j 5 Qyh^i}) u ■ (8) 

Here the index U denotes an expectation value in the path integral over gauge variables. A 
similar formula is also valid for scalar mesons, where 75 is replaced by the unit Dirac matrix. 
Since in QCD the correlation analogous to (§), apart from the factor 2 in front of the connected 
contribution, belongs to rj and 77' mesons, one can call the corresponding boson a-eta. 

As discussed in ref. @, the expectation value in eq. @ can also be expressed by a "noisy 
estimator" in terms of the n bosonic pseudofermion fields ((f)). In the present case we have 



4 £ + (<t>*M] [(<vv<) + - 2^(0*^,)) , (9) 

j,k=i 



with the notation 



>jx = J2(^Q + f'.Or./.'j,, ■ (10) 



In the last formula /ij (j = 1, . . . , n) denote the real parts of the roots of the polynomial used 
in the local bosonic updating algorithm. 

The advantage of the noisy estimator in eq. (0) is that in principle it can be obtained with 
very little work from the pseudoscalar fields during unquenched updating. Unfortunately, this 
does not work in practice: the result is too noisy and the signal cannot be obtained in this way 
0. What remains is the evaluation of the fermion propagators appearing in eq. (^. 

The "connected" contribution in the second term has been discussed in the previous subsec- 
tion. The behaviour of the first "disconnected" contribution on 4 3 • 8 lattice at /3 = 2.3, K = 
0.180 is shown in figure [5[ Compared to the connected term, the disconnected contribution is 
two to three orders of magnitude smaller. On the other hand, the disconnected contribution 
decreases considerably slower as a function of the distance. It shows an effective mass of the 
order of m^ sc ~ 0.8, whereas the effective mass in the connected contribution (the a-pion mass) 



is about m conn ~ 1.3 . This behaviour is qualitatively similar to the one observed in QCD [13 



As argued in ref. [IB], in the quenched approximation one can also extract a "quenched" 



a-eta mass from the ratio of the disconnected to connected contributions. Denoting the ratio 
of the two terms at timeslice distance t by R(t), one has 

2 

TfL 

R(t) = const. + t — 2_ . (11) 
2m,,- 

Here m v and m n stand generically for the masses of the flavour singlet and flavour non-singlet 
pseudoscalar mesons, respectively. In our case we get from (|TT| ) a quenched a-eta mass of 
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about m a _ eta ~ 0.1, although the linear fit on our small 4 3 • 8 lattice is not very good. Of 
course, we are finally interested in numerical results with dynamical gluinos, where the linear 
behaviour in ([11]) is irrelevant, and the a-eta mass can be obtained from a usual kind of fit of 
the combination in eq. (§) according to eq. (0). The results shown in figure |5| were obtained 
by direct evaluation of the necessary inverses on 100 independent configurations. In fact, for a 
more precise calculation of the disconnected contributions, especially on larger lattices, there 
are better methods known (see ref. [PH ). 

In summary: the quenched calculations give useful information on the spectral properties of 
fermion matrices and the range of hopping parameters relevant for the continuum limit, which 
help in setting up numerical simulations with dynamical gluinos. 
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Figure 1: The quenched distribution of the lower spectrum boundary p min of 
Q*Q on 8 3 • 16 lattice at f3 = 2.30, K = 0.20 . 
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Figure 2: The quenched distribution of the lower spectrum boundary \R m i n 
of Q R = (Q + Qt)/2 on 8 3 • 16 lattice at (3 = 2.30, K = 0.20 . 
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lattice: 8®8®8®16 
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Figure 3: The a-pion mass squared as a function of the inverse hopping 
parameter on 8 3 • 16 lattice at (3 — 2.3 . The line is a linear fit used to determine 
K 
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lattice: 8<8>8®8®16 
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TT a and cr a mass (open = 4 3 <S>8) 
13 = 2.3 



[] 
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1/K 

Figure 4: The a-pion and a-sigma mass as a function of the inverse hopping 
parameter on 8 3 • 16 lattice at f3 — 2.3 . The open symbols refer to 4 3 • 8. The 
curves show linear fits for the mass squares. 
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lattice: 4040408 
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Figure 5: The disconnected a-eta correlation as a function of the timeslice 
distance on 4 3 • 8 lattice at (3 = 2.3, K = 0.180 . The curve is a fit according; 
to (0) for < t < 4, giving a mass mdi SC = 0.84(15). 
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